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Многочисленные теоретические и экспериментальные работы [4, 5, 
6] посвящены исследованию зависимости скорости обмена воды и нефти 
в единичном блоке породы. На перемещение фронта вытеснения могут 
оказывать влияние как свойства матрицы и трещин, так и физические 
свойства жидкостей. В случаях очень высокой вертикальной 
проницаемости трещин фронт вытеснения в упрощенной модели 
принимается горизонтальным. 

Различные авторы используют в своих экспериментальных работах 
модели, характеризующиеся различной геометрией блоков, формой и 
свойствами контактирующих с трещинами поверхностей, разными 
скоростями перемещения уровня воды в трещинах (эквивалент скорости 
погружения блока), а также разной интенсивностью проявления сил, 
участвующих в процессе вытеснения, и т. д. 

Выводы. Во всех экспериментах по вытеснению как при 
преобладании гравитационных или капиллярных сил, так и в тех случаях, 
когда ‘эти силы вместе участвуют в процессе вытеснения, кривая 
зависимости нефтеотдачи от времени подобна кривой, приведенной на 


рис. 1, 6. 


Типы кривых зависимости нефтеотдачи К от времени: - 
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Кривая на рис. 1, 6 допускает некоторую линейную аппроксимацию 
в полулогарифмических координатах для промежуточных значений 
времени. 

Основываясь на теоретическом рассмотрении и результатах 
экспериментальных исследований Маттакса, приведенных в безразмерном 
виде, можно прийти к следующим выводам. 


а. При малой высоте блоков, когда капиллярная пропитка определяет 


а 
процесс вытеснения, время вытеснения пропорционально квадрату высоты 
а блока Н 
а 6. В высоких блоках, где гравитационные силы управляют 


процессом вытеснения, время вытеснения пропорционально высоте 
блока Н. 
в. В случае «а» нефтеотдача может быть представлена как функция 
15, а в случае «б»—как функция 5. 
Используя зависимость нефтеотдачи от логарифма времени, можно 
сформулировать положение о критериальном времени, необходимом для 
| получения 50% - ной отдачи, — 15. Сопоставление величин {0 ѕв, 
полученных в экспериментах, проведенных при различных условиях 


вытеснения, может оказаться полезным при прогнозировании и анализе 


‚ "ж: ~ разработки реальных залежей. 


Критерии подобия процессов несмешивающегося вытеснения в 


МГ пористых средах были детально рассмотрены Раппопортом [8]. Основные 


соотношения, помогающие переносить результаты лабораторных 


экспериментов на натурные условия, позволяют предсказывать поведение 


трещиноватых пород-коллекторов. Учитывая, что градиенты давления в 


трещинах незначительны по сравнению С гравитационными И 


капиллярными силами, необходимо, чтобы критерии подобия учитывали 
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только эти силы. Если найдены правильные критерии подобия, то 
использование в экспериментах образцов, отобранных из данной залежи, 
позволяет воспроизводить процессы, протекающие в этой залежи. 

Критерии подобия надежны только в том случае, когда соблюдаются 
основные допущения относительно однородности матрично-трещинной 
системы. Однако при существенной экстраполяции результатов 

„ экспериментальных работ следует соблюдать целый ряд 
предосторожностей. 

И Требования теории подобия. Если при помощи модели 

предпринимается попытка предсказать поведение залежи, образованной 


‚блоками с известными характеристиками, необходимо удовлетворить 


некоторому ряду основных требований. Эти требования основаны на 


классических принципах, развитых Раппопортом [8] для процессов 
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несмешивающегося вытеснения. А 


Вопросы, связанных с экспериментами по капиллярной пропитке 


Ф 


Основными факторами, влияющими на результаты экспериментов и 
2 _ ограничивающими их использование, являются неоднородность и тип 
А пород, тип и свойства жидкостей, а также условия проведения . 


пе Р", эксперимента. 


а. Неоднородность пород 
В пределах залежи неоднородность пород ограничивает 
применимость зависимости нефтеотдачи от времени, получаемой из 
экспериментов Маттакса [4] (рис. 2). Парсон [10] изучал влияние : 
а . неоднородности в опытах по капиллярной пропитке на образцах размером 
5,08х 5,08 х 2,54 см (рис. 2, а), также исследуя зависимость нефтеотдачи от Р 
я ° времени. Затем образец разрезали на четыре кубика, и опыты проводились 


на каждом из них. Результаты, полученные для образца /, дали заметный 


е 


разброс значений нефтеотдачи (рис. 2, 6), в то время как результаты 


опытов на кусочках из образца 2 указывали на значительную однородность 


породы (рис. 2, в). 


Влияние неоднородности коллектора на зависимость 


нефтеотдачи от ѓ;.„ [10] 
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а – форма образцов; 


‚ _ результаты экспериментов: б – для неоднородных пород, в – для 
однородных пород. / - размеры образца 5,08 х 5,08 х 2,54 см; 
2 – 5 – кубики соответственно № 1—№ 4 с ребром 2,54 см 


Рис. 2 


6. Литологические свойства 


Опыты, проведенные при одинаковых условиях, могут дать разные 


результаты из-за различия минералогических характеристик пород [11]. 


Результаты экспериментов на различных литологических разностях (типы 


А, и Аз) схематически показаны на рис. 3, а. При разных значениях С С А 


расхождение результатов может целиком и полностью быть прямым 


следствием вариаций литологического состава (рис. 3, 6). 


Наблюдается существенное различие зависимостей нефтеотдачи от 


времени для разных литологических разностей при самопроизвольной 
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капиллярной пропитке (см. рис. 3, 6). Наблюдаемую тенденцию можно 


объяснить тем, что возрастающая карбонатность пород определяет их 
более высокую гидрофобность и соответственно меньшую нефтеотдачу. 
в. Природа жидкостей 


Эксперименты, проведенные при одних и тех же условиях на 


ат 
подобных друг другу образцах, но с применением различных по свойствам 


жидкостей, показали, что кривые зависимости нефтеотдачи от времени 


ПИ 


отметить значительные расхождения между опытами с обычной морской 


существенно расходятся Результаты экспериментов позволяют 


"водой и морской водой с добавлением ингибиторов. При использовании 
искусственной морской воды получены непредсказуемые результаты. 
Резкое снижение нефтеотдачи в опытах с искусственной морской водой 
объясняется очень быстрым снижением межфазного натяжения на ее 
границе. Это свидетельствует о том, что в случае изменения межфазного 
натяжения с течением времени результаты лабораторных экспериментов 
становятся полностью непригодными для моделирования реальных 


процессов 


г. Влияние положения поверхности блока, контактирующего 


с вытесняющейся жидкостью 


Представляется необходимым сравнить между собой результаты 
процесса вытеснения при разном положении граней, поскольку это влияет 
на тип и направление процесса вытеснения. Достаточно большое число 
представительных случаев было исследовано Иффли (см рис. 4). Образцы, 
боковые поверхности которых были герметизированы, в случаях 4, ОП и Е 
имели открытые верхний и нижний торцы, а в случаях В и С — только 


один торец. Несмотря на то что величина С С А поддерживалась во всех 


опытах на одинаковом уровне С С К = 150, во всех пяти исследованных 
вариантах изменение нефтеотдачи во времени существенно отличалось из- 


за различного положения проницаемых поверхностей. 


Влияние литологического состава пород на зависимость нефтеотдачи 


от времени #5. [11] 
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а – типы капиллярной пропитки; б - зависимость нефтеотдачи от 


времени для различных пород; в — минеральный состав образцов 


Рис. 3 


Из результатов экспериментов, представленных на рис. 4, следует: 
– в эксперименте В реализованы наиболее неблагоприятные условия, 
так как нефть вытеснялась противотоком и в направлении, обратном 


действию гравитационных сил; 


Зависимость нефтеотдачи от безразмерного времени #5. 
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пять различных положений плоскости ВНК при одном и том же 
значении капиллярно-гравитационного отношения С С А = 150 [11] 


Рис. 4 


— в эксперименте С реализованы более благоприятные условия для 
вытеснения нефти, так как несмотря на противоточную капиллярную 
пропитку направление вытеснения совпадало с направлением действия 
гравитационных сил; 

— наиболее благоприятные условия реализованы в случае Е, когда 
все силы (гравитационные и капиллярные) помогали вытеснению нефти; 

— лучшие результаты вытеснения получены в случае А (по 
сравнению с 20) вследствие меньшей эффективности противоточного 
вытеснения даже при благоприятном направлении действия 
гравитационных сил по сравнению с прямоточно-противоточным 


капиллярным вытеснением (случай А). 


д. Влияние вязкости 


Влияние вязкости нефти на процесс вытеснения изучено на 
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численных математических моделях [12] для пяти различных случаев (Ин 


= 0,01; 1; 5; 30; 100 мПа с) при постоянной вязкости воды, равной 1 мПа с. 


При противоточной капиллярной пропитке распределения 


насыщенности существенно отличались при различных вязкостях нефти. 


Причем было показано, что только при низких вязкостях нефтей 


реализуется фронт вытеснения. При высоких значениях вязкости 


выраженный фронт вытеснения отсутствует, а наблюдается постепенное 


увеличение водонасыщенности по длине модели (рис. 9), 


7 Влияние вязкости нефти на скорость капиллярной пропитки 
отражено на рис. 6 для величин нефтеотдачи, больших 6%, причем 
зависимость нефтеотдачи от времени приближенно описывается прямой 
линией в логарифмических координатах. При более низкой вязкости нефти 
за одно и то же время достигается более высокий коэффициент 
нефтеотдачи. Иное отношение приведено на рис. 7, где иллюстрируется 
зависимость величины относительного времени нефтеотдачи как функция 
ВЯЗКОСТИ. 

Под относительным временем нефтеотдачи понимается отношение 


(1) нефтеотдачи на данный момент времени при данном значении Ин 
к нефтеотдаче к этому же моменту времени при ди = 1 мПа с, которое 


показывает, насколько увеличивается время, необходимое для извлечения 
того же количества нефти при большей ее вязкости. Расчеты проводились 
для образца с площадью сечения А =7,92 см”, длиной [ = 50 см, 
пористостью Ф = 32%, начальная водонасыщенность принималась равной 


9.2% при неизменной вязкости воды = 1 мПа с. Ра 
р В е 


Ж». 
в: 


Распределения насыщенности в различные моменты времени 


для нефтей разной вязкости [12] 
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5в.НАЧ = 9,2 %; [= 50 см; А = 7,92 см; Ф = 32,1 %; дв = 1 мПа с; К = 0,2 мкм? 
Рис. 5 


е. Влияние характеристик пород 


Капиллярное давление зависит от изменения пористости, 
проницаемости и распределения пор по размерам. Из результатов 
наблюдений известно, что относительная проницаемость и начальная 
водонасыщенность играют второстепенную роль в процессах пропитки, 


так что капиллярные силы при капиллярной пропитке в основном 


являются функцией величины ү К/Ф. 


Зависимость нефтеотдачи от времени капиллярной пропитки для 


нефтей различной вязкости 
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Рис. 6 


Изменение безразмерного времени (‚т в зависимости от вязкости 


нефти 
быт 
ю0 
ю |- 
7 
01 Ч 
и п 1 10 10 


Вязкость нерти, м Па-с 


ыт = отношение времени вытеснения заданного объема нефти некоторой ВЯЗКОСТИ 


к времени вытеснения того же объема нефти вязкостью 1 мПа с [12] 


Рис. 7 
ж. Влияние подвижности М =л.> 


В экспериментах Дю-Прея [3] было установлено, что при изменении 


подвижности М = л, от 0,025 до 0,9 зависимость нефтеотдачи от 


® 


безразмерного времени несколько изменяется. 


з. Влияние капиллярно-гравитационного отношения СС К = л; 


Капиллярно-гравитационное отношение показывает, какая из двух 


сил преобладает в механизме вытеснения, т. е. какая сила оказывает 


наибольшее влияние на вид зависимости нефтеотдачи от времени. 

= Эксперименты были проведены Иффли [11] на образцах высотой от 
5 до 200 см и проницаемостью от 10:10° до 1 мкм? с применением 
реальных флюидов и пород примерно одного и того же минерального 
состава (силикаты = 60%, глинистые = 15%, карбонаты = 25%). Вид 
‘эксперимента (обозначаемый индексом А.1) соответствовал прямоточно- 
противоточному вытеснению нефти через верхний торец модели. 
Изменяющаяся величина С С А: 

п -сав-2/8УК/Ф МК/Ф (1) 
Дрећ ео 


составляла от 20 до 1000. Результаты экспериментов, приведенные на 
рис. 8, показывают со всей очевидностью, что нефтеотдача строго 
определяется величиной фактора С С А независимо от других условий 
эксперимента. Основной вывод из этих опытов заключается в том, что 


вытеснение нефти происходит тем быстрее, чем ниже значение С С К, что 


в свою очередь указывает на то, что возрастание гравитационных сил 


ч = 
А 


ведет к ускорению вытеснения. 


Аналогичные эксперименты были выполнены на образцах, имевших 


форму параллелепипеда с одной проницаемой боковой поверхностью [9]. 
Результаты, представленные в виде зависимостей нефтеотдачи от /5 и (5, 


имеющих одинаковый вид, приведены на рис. 8. Из них следует, что при 


меньших значениях ССА что отвечает возрастанию роли 


е 


гравитационных сил в процессе вытеснения, отмечается более быстрое 


извлечение нефти. 


Зависимость нефтеотдачи от ѓ„ для различных значений капиллярно- 


гравитационного отношения С С К [11] 
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Рис. 8 


3. Процессы гравитационно-капиллярного вытеснения 


при их моделировании на центрифуге 


Эксперименты Кайта [7] е учетом критериев подобия по оценке 


нефтеотдачи при действии гравитационных сил представляют собой 


быстрый и относительно точный метод. Поэтому их широко используют 
при подготовке данных для проектирования разработки. 

Как упоминалось ранее, если в процессах капиллярной пропитки для 
учета влияния капиллярных сил использовать численный метод Блэйра 
[12], необходимо точно знать зависимости капиллярного давления и 
относительных проницаемостей от насыщенности. Метод 
центрифугирования не требует знания этих функций, во всяком случае, в 
тех экспериментах, в которых используются представительные образцы. 

Кайт замечает, что функции относительных проницаемостей очень 
часто экспериментально определяются при высоких скоростях 
фильтрации, при которых капиллярные и гравитационные градиенты 
давления пренебрежимо малы по сравнению с градиентами фильтрации. 
И если данные таких экспериментов затем используются в численных 
расчетах работы реальной залежи, то полученные результаты будут сильно 
отличаться от действительных показателей разработки. 

Помимо этого положительной стороной экспериментов на 
центрифуге является то, что концевые эффекты могут быть значительно 
лучше смоделированы на центрифуге, чем при использовании численных 
моделей, в которых в качестве граничных насыщенностей блоков могут 


быть приняты сомнительные величины. 


=-У 


а. Влияние капиллярно-гравитационного отношения СС К= л; 
А . Ж 


м. ; 


Значение гравитационных сил было выявлено в экспериментах 
Маттакса, в которых было показано, что для блоков матрицы высотой 1,82 
м гравитационные силы оказывают в 12 раз более эффективное влияние на 
процесс вытеснения, чем капиллярные. 

В экспериментах, проведенных Дю-Преем [9], сравнивались 


обычная капиллярная пропитка и центрифугирование для различных 
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величин С С А. 

Непосредственное влияние на вид зависимости нефтеотдачи от 
времени оказывает величина ССА. Причем анализ результатов 
показывает, что их можно упорядочить с учетом безразмерного 


"капиллярного" времени &%, или безразмерного "гравитационного" 


времени бгр. 
> 7 . По отношению к безразмерному времени {5 к одна и та же величина 
нефтеотдачи достигается быстрее при снижении С С К, т. е. в тех случаях, 
когда роль гравитационных сил в процессе вытеснения увинчивается 


(рис. 9, а), но при этом результирующие кривые для разных значений 


ССР расположены дисперсно. 

На зависимость нефтеотдачи от ССА фактор {5 оказывает 
аналогичное влияние: одна н та же величина нефтеотдачи достигается 
быстрее при возрастании гравитационных сил, но результирующие кривые 
располагаются значительно ближе одна к другой (рис. 9, 6). 

В силу дисперсии кривых "нефтеотдача — время" 
в билогарифмических координатах Дю-Прей предложил анализировать 
результаты экспериментов как функцию среднего времени вытеснения. 

г 7 ә. Среднее время вытеснения равно безразмерному времени, при котором | 
нефтеотдача составляет 50% (см. рис. 9). Это время обозначается {со 5р. 
ан" Изучение зависимости С С А при вариациях (50.5р в интервале от | до 
10 показало, что в билогарифмических координатах между ними 
существует линейная связь. Причем наклон прямых С С К зависит от того, 
изображаются ли они как функции {к или как функции (9.55. „ | 
Из рис. 10 видно, что линия зависимости СС К от {созк идет круче 
в случае преобладания гравитационных сил, чем в случае преобладания 


капиллярных. 


Зависимость нефтеотдачи от безразмерного времени %,, (а) и (б) [9] 
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1 ~ капиллярная пропитка; 2 ~ гравитационное вытеснение [9] 


Рис 10 
б. Влияние геометрических размеров блока (высоты и диаметра) 


Влияние различной высоты блоков на среднее время вытеснения 
после получения результатов лабораторных исследований с учетом 
критериев подобия было перенесено в натурные условия. Полученные из 
зависимости нефтеотдачи от безразмерного времени размерные величины 
среднего времени вытеснения (55) были исследованы как функции 
высоты блоков (рис. 11) с учетом того, что эксперименты Дю — Прея 
ВЫЯВИЛИ линейную зависимость исследуемых величин 


в билогарифмических координатах. Нарушение линейной зависимости 


наблюдается для блоков, высота которых превосходит 1 м. Это 
соответствует тому, что для блоков высотой меньше 1 м, 15 
пропорционально Н Прямая пропорциональность при высотах блоков 
Н>1 м может быть нарушена в результате методических трудностей при 
проведении экспериментов на центрифуге при высоких значениях 


центробежного ускорения а > 10 г. 


Изменение среднего времени вытеснения / 5к в зависимости от 


высоты блоков Н, наблюдавшееся в экспериментах Дю-Прея [9] 


Нем 
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Рис. 11 


В экспериментальных работах Дю – Прея при постоянстве С С А 
было обнаружено более быстрое вытеснение нефти из блоков большого 
диаметра. Но результаты этих экспериментов не могут рассматриваться 
в качестве окончательных, и для выводов требуются дальнейшие 


исследования. 


Резюмируя все выше сказанное, можно сделать следующие выводы: 


а. Получена удовлетворительная корреляция между теоретическим 


тае У 


анализом и экспериментальными исследованиями процессов капиллярной 
пропитки. 

6. Процессы вытеснения в малых блоках в основном определяются 
капиллярными силами, а в больших (высоких) блоках — гравитационными 
силами. : 

в. Если в процессе вытеснения преобладают капиллярные силы, то 
время вытеснения пропорционально квадратному корню из высоты блока, 
при преобладании же гравитационных сил время вытеснения 
пропорционально высоте блока. 

г. Зависимость нефтеотдачи от безразмерного времени к и /5, 
полученная в экспериментах, так же оказалась очень близкой 
к результатам теоретических исследований. Естественно, что результаты 
этих исследований оказываются менее коррелируемыми при 
использовании неоднородных по литологическим свойствам образцов и 
в том случае, когда свойства жидкостей оказывают влияние на результаты 
экспериментов при взаимодействии с пористой средой модели. 

| д. Моделирование блоков большой высоты путем использования 
метода центрифугирования оказалось очень быстрым и надежным 
методом, но, однако, при высоких частотах вращения ротора надежность 
результатов моделирования неудовлетворительна. 

е. Значительный шаг в направлении надежного прогнозирования 
нефтеотдачи из единичного блока матрицы удалось сделать благодаря 
результатам, полученным на основе моделирования капиллярной пропитки 
как обычным методом, так и методом центрифугирования. в 

ж. Критическое отношение к физическому моделированию и замена 
его численными методами на первом этапе казалась логичной. Однако 
невозможность учета в численных методах реального вида функций 


относительных проницаемостей привела к тому, что в очень редких 


БА 


О?У 


ав 
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случаях удавалось получать надежные результаты при одновременном 
проявлении гравитационных и капиллярных сил. | 

з. Сложность явлений, которые составляют процессы капиллярной 
пропитки, заставляет предполагать, что использование обеих процедур: 
экспериментальной (с учетом критериев подобия) и численной – в 


дальнейшем может привести к созданию более общей теории. 


4. Другие аспекты механизма вытеснения нефти из единичного блока 


Дренирование нефти из блока за счет сил гравитации наблюдается 


в тех случаях, когда расширяется газовая шапка. Характеристики процесса 


` дренирования единичного блока можно рассчитать, когда известны 


размеры блока и свойства жидкостей. 

Очень важным представляется сопоставление высоты блока Н; 
с высотой капиллярного подъема Н, и высотой прорыва А, которые 
определяются по кривой капиллярного давления. 

В низкопроницаемых блоках небольшой высоты в силу того, что 
высота капиллярного поднятия больше высоты самого блока, вся нефть 
должна быть удержана в блоке. Тогда же, когда высота блока больше, чем 
высота капиллярного поднятия, конечная нефтеотдача будет получена 
к моменту наступления равновесия между гравитационными 
и капиллярными силами. При этом нефть вытесняется из зоны, 
заключенной между нижней гранью блока и уровнем, соответствующим 


высоте капиллярного поднятия, т.е. из зоны капиллярного поднятия. 


Следовательно, только из блоков очень высокой проницаемости, 


характеризующихся низкими капиллярными давлениями, или из очень 
высоких блоков может быть получено существенное количество нефти. 
Таким образом, высота блока и характерное капиллярное давление — 


наиболее важные параметры с точки зрения эффективности разработки. 


4 
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Численное моделирование процесса выполнено Саиди [16] для 
блоков высотой Н. = 3 м, Др = 0,6 г/см, ин = 1 мПас, К(Н)5видч 
1,26-10° мкм”, Ф(1- види) = 0,083 и Кы= (5,)' [7]. 

Нефтеотдача очень быстро возрастает до Ко, особенно для низких 


значений конечной нефтеотдачи (сравните кривые для Ах = 30% 


и Коо = 50%0 от объема извлекаемой нефти). Результаты моделирования 


приведены на рис. 12. 


Зависимость нефтеотдачи от времени при вытеснении нефти газом 


газовой шапки [16] 
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Рис. 12 Е 


_ Гравитационное дренирование можно изучать с использованием 
методов центрифугирования, при этом целый ряд условий должен быть 
принят во внимание. Так, максимальный гравитационный потенциал будет 
реализован в верхней, а не в нижней части образца, а это требует 
использования специальных кернодержателей, позволяющих исследовать 
образцы, окруженные газом. 

| Добыча при нагнетании воды. Вопросы, связанные с нагнетанием 


воды, изучались целым рядом исследователей [9, 13, 14, 15]. И во всех 
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работах, как теоретических, так и экспериментальных, ставилась цель 
оценить зависимость нефтеотдачи от времени и доли воды в потоке 


в зависимости от накопленного объема закачки. Исследовалось влияние на 


процесс различных параметров, таких как отношение проницаемостей 


трещин и матрицы, действующих градиентов давления и т.п. 


ре -. 5. Теория нефтеизвлечения из единичного блока 


МЕ В предыдущих разделах рассматривались процессы извлечения 
- - нефти из единичных блоков, причем особое внимание обращалось на 
физические аспекты вытеснения нефтей при капиллярной пропитке, т.е. 


`при вытеснении нефти водой, или при дренировании, т.е. при вытеснении 


нефти газом. 


Единичный блок при этом рассматривался как однородный, 

и поэтому весь процесс вытеснения описывался уравнениями, пригодными 

для описания процессов фильтрации в однородных породах. 

С В действительности задача значительно более сложная, так как 
в трещиноватых породах-коллекторах блоки матрицы могут обмениваться 
жидкостями только через прилегающие к ним трещины. Таким образом, 

К вдоль направления вытеснения вытесняющий агент стремится 

` перемещаться по трещинам с большей скоростью, чем по блокам, 

благодаря высокой неоднородности, создаваемой наличием сети трещин. 
С другой стороны, непрерывное проникновение вытесняющей жидкости из 
трещины в блоки приводит к снижению разницы в скоростях движения 
вытесняющего агента в обеих средах: трещинах и пористых блоках. 
Взаимная обусловленность процессов вытеснения нефти из блоков 
и трещин требует нового подхода к трактовке вопросов, связанных с 


вытеснением нефти из блоков матрицы в трещиноватых породах- 


коллекторах. Этот подход основан на рассмотрении трещиноватых пород- 


е 
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коллекторов как целого. Так, аналитическое исследование трещиноватых 
8 | пород-коллекторов основано на континуальном подходе ("вложении" друг 
в друга сред с различной пустотностью) в формулировке Баренблатта [16]. 
|. В численных моделях, тем не менее, возможно представление 
трещиноватого коллектора в дискретном виде, т.е. в виде отдельных 


трещин и блоков матрицы. 
ањ < 


4 


5.1. Континуальное приближение Баренблатта 


2. 


2009 Континуальное приближение Баренблатта [16] для случая 
однофазных течений в трещиноватых породах-коллекторах обобщено на 
‘случай двухфазных течений несмешивающихся жидкостей. Уравнения 
фильтрации выписываются для течения каждой фазы через составную 
среду трещиноватого резервуара — трещины и пористые блоки; переток 
вытесняющей жидкости из трещин в блоки и истечение вытесняемой 


жидкости из блоков в трещины описывается членами источник — сток 


44 Е 
в уравнениях сохранения массы в соответствии с законом, определяющим 
вытеснение флюида из блока, т.е. функцией капиллярной пропитки. 
а Боксерман, Желтов и Кочешков [6] представили уравнения 
| ~. < Сохранения массы воды и нефти в трещинах следующим образом: 
25 05 ее 
Е | Фу 22 + дуи) +ир = 0; (2) 
а $ - КЕС | 01 
08р ` . ҹ : 
НЭ: Е С 
Ф! + дїуирә +ин = 0, + (3) 
дї у 
где Ови Он - соответственно скорости воды и нефти из уравнений Дарси: : 
р 
4 ЕЕ; 
| — _ К, Ков2(5в2) 
22150-82082 
ир = — Рғаау; 
Ив г, 
Е т < К, Кон2(5н2) 
2 опорно НО $ 
ана = ЕАН раду, (4) 
г Е Д Ан 5 
Ей е7 Ы ить. 
& Е Ў У ‹ 
: > 
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Поскольку объем трещин очень мал, допускается, что определенный 


объем воды, входящей в трещины, проникает вследствие капиллярной 


пропитки в блоки. При допущении несжимаемости жидкости объем воды, 


впитываемой блоками, равен объему извлекаемой из трещин нефти и, 


* 


ж ж 
следовательно, ин =-ин = |и |. 


‚ * Это предполагает определение необходимых функций на основе 
.. экспериментальных исследований. Определив границу поверхности, 
- ограничивающей объем /(ї*) впитываемой воды, через 5 (х, х›, хз, И 


можно получить значения искомой функции {*(х,, х›, хз) из следующего 


интегрального уравнения, описывающего баланс воды в трещиноватом 


коллекторе: 
, | | 25) А = ад) | (5) 


Функция переноса, предложенная Боксерманом и др. [6], имеет 


ё следующий вид: 
є. АІ 
7 А 5ссоѕө 4 Ку/Ф, [ 52 соз0/К/Ф | 2 
и (0) = Ф 5 — 1 , (6) 
2 Ан Ин 
Уч ‘>. = 
1" где А — постоянный коэффициент. 
ей 
5. Боксерман с сотрудниками [6] дал решение для одномерного 
вытеснения нефти водой. В этом случае задача сводится к решению 
И: 
интегрального и дифференциального уравнений: { 
И Јас 0) саи 
и 0 $ 
7. ї 053 Ф ссоѕ0 ү Ку/Фу 0$, » * 
15-40) (5) 8 20|, -г(,)]е0, 0) 
Е х Ин 16 
к = > ечи ~ в 
Е + 9 в: & 
4 1 | , А 
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где: 
Х=Х/Ё 
1; =10с0504[Ку1Ф 52° / ин; 
Е(5в)= Ко.в(5в)/[Ко.в(5в)+ (ив / ин Ко н(5в}}; 
и" (#;)= (А/2)Ф, Ку; [6 соз0 5? [857% 52 595 /ин| 


(8) 


где Г – характерная длина, например средний размер блока. 
Пропитка определенного объема блока происходит за время Е, равное 
15, удовлетворяющее следующему условию: 


1—2 ()<2,. 


Решение уравнений (7) при постоянном расходе д = сопѕі 


и начальных условиях Х; = 0 при &=0 имеет вид: 


бах, (9) 
где: 
2 
МЕ тА 52ссоѕ0 |К, 1Ф, Ф, 5р (10) 
44 Ин | 


5. Теоретический подход Брестера 


Брестер [15] рассматривал уравнения совместного течения воды 
и нефти в трещиноватых породах — коллекторах как функцию источника, 
описывающую вытеснение нефти из блоков матрицы за счет процессов 
капиллярной пропитки, гравитационного перераспределения и действия 
градиента давления в трещинах. Кроме того, модель допускала течение 
жидкостей через блоки без вытеснения, т. е. обмен одними и теми же 
жидкостями между блоками и трещинами, так же как в случае 


однофазного течения. Таким образом, течение возможно только на тех 
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участках блоков, где величины насыщенности допускают наличие 
подвижных флюидов. Исходя из этой точки зрения, течение какой-либо из 
жидкостей в трещинах и блоках можно рассматривать как течение 
в модели, представленной системой трещин, параллельной системам 
блоков, насыщенным рассматриваемыми флюидами. Возрастание 
насыщенности блоков каким-либо из флюидов аналогично увеличению 
| эффективной проницаемости трещин, т.е. эффективные проницаемости 
для каждой из жидкостей в трещиноватой системе эквивалентны средней 
проницаемости трещиноватого коллектора (системы трещин и блоков) 
и относительной проницаемости, зависящей от насыщенности пустот 
обоих типов: блоков и трещин. 


Для системы параллельных трещин и блоков средняя проницаемость 


среды составит 
К =(К, В, +К, В,)/(В, + В,), 


где К. и В соответствуют проницаемости и толщинам рассматриваемых 


сред. Так как 


в некоторых коллекторах, характеризующихся блоками большой 
протяженности, вклад блоков (К›В›В,) в суммарную проницаемость 
весьма значителен. Тем не менее, для коллекторов, характеризующихся не 
пренебрежимо малыми проницаемостями блоков матрицы, также пригодна 


приведенная модель. 
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